Edge Dynamics in a Quantum Spin Hall State: 
Effects from Rashba Spin-Orbit Interaction 



O 

(N 

(N 



CD 



o 

> 

\^ 

^' 

O 

o 



X 



Anders Strom^, Henrik Johanncsson^ , and G. I. Japaridze^'^ 
^ Department of Physics, University of Gothenburg, SE 4.12 96 Gothenburg, Sweden 
^ Andronikashvili Institute of Physics, Tamarashvili 6, 0177 Tbilisi, Georgia and 
Ilia State University, Cholokashvili Avenue 3-5, 0162 Tbilisi, Georgia 

We analyze the dynamics of the hehcal edge modes of a quantum spin Hall state in the presence 
of a spatially nonuniform Rashba spin-orbit (SO) interaction. A randomly fluctuating Rashba SO 
coupling is found to open a scattering channel which causes localization of the edge modes for a 
weakly screened electron-electron (e-e) interaction. A periodic modulation of the SO coupling, with 
a wave number commensurate with the Fermi momentum, makes the edge insulating already at 
intermediate strengths of the e-e interaction. We discuss implications for experiments on edge state 
transport in a HgTe quantum well. 
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The quantum spin Hall (QSH) state is a recently dis- 
covered phase of electronic matter in two dimensions 
which exhibits no symmetry breaking P, Q . Being an 
example of a topological insulator the QSH state is 
gapped in the bulk but exhibits massless edge modes. 
These modes propagate in opposite directions for oppo- 
site spins, each helical mode having a time-reversed copy. 
In the range of electron-electron (e-e) interactions probed 
in experiments elastic backscattering induced by a 
single (nonmagnetic) impurity or by quenched disorder 
is suppressed by time- reversal invariance 0, Q . For sam- 
ples of a length smaller than the spin dephasing length Q 
one thus expects that the edge modes are protected, re- 
sulting in nondissipative transport with a quantized Lan- 
dauer conductance jh. This expectation is supported 
by experiments holding promise for exploiting QSH 
edges in future integrated circuit technology as a means 
to reduce power dissipation as devices become smaller. 

The question whether the QSH edge is stable also 
against e-e interactions in the presence of a Rashba spin- 
orbit (SO) interaction has received less attention. One 
should here recall that SO interactions — which allow for 
spin-flip scattering while respecting time-reversal invari- 
ance — come in different guises. In the case of a HgTe 
quantum well — the best studied platform for realizing 
a QSH state 0, B ~ the intrinsic SO interaction in the 
atomic p orbitals crucially contributes to the very special 
properties of this material: For a quantum well (QW) 
thicker than a critical value it inverts the band struc- 
ture, pushing the p band above the s band. A simple 
model calculation reveals the presence of a single pair 
of helical edge bands inside the inverted bulk band gap, 
showing that this intrinsic SO interaction is what makes 
possible the QSH state This is to be contrasted 

with the Rashba SO interaction which originates from the 
gate-controllable inversion asymmetry of a QW • Its 
strength can be huge in the bulk metallic phase of an in- 
verted HgTe structure, with an effective coupling several 
times larger than for most other known semiconductor 



heterostructures In a QSH state with an insulating 
bulk, charge and spin are transported by the edge modes. 
With the bands for these modes dispersing all the way 
from the bulk valence band to the bulk conduction band 
0, one expects a sizable Rashba interaction to still op- 
erate at the edge of the sample, making it important to 
inquire about its effect on the edge modes. 

In this Letter we analyze the dynamics at the QSH 
edge in the presence of a Rashba SO interaction. Treat- 
ing the Rashba interaction as a perturbation, we find that 
it has no effect on the propagation of edge modes when 
it is spatially uniform. The more realistic case when the 
Rashba coupling fluctuates randomly in space 1^ |13[ is 
different, however. Together with a weakly screened e-e 
interaction, the Rashba interaction now opens a scatter- 
ing channel for the edge modes which may cause them to 
localize. For a sufficiently strong e-e interaction we find 
that the localization length is shorter than the length of a 
typical sample Q , implying that the edge modes do local- 
ize, causing the QSH state to collapse. As the screening 
in a QW can be controlled by varying the thickness of the 
insulating layer between the well and the nearest metallic 
gate, our prediction is amenable to a direct experimental 
test. Interestingly, by subjecting the edge to a periodi- 
cally modulated Rashba interaction, with a wave num- 
ber commensurate with the Fermi momentum, we find 
that the edge becomes insulating already at intermediate 
strengths of the e-e interaction. A periodic and tunable 
Rashba modulation, produced by a sequence of equally 
spaced nanosized electrical gates, may thus be used as a 
current switch in a QSH-based spin transistor. 

To model the edge dynamics of a QSH state we use 
two electron fields ^^{x) = ip^{x)e^''''^ and Vl'^(x) = 



where l annihilates a clockwise 



[counterclockwise] propagating edge electron with spin- 
up [spin- down], quantized along the growth direction z 
of the QW. The fields if}^ and t/j^ are slowly varying, with 
the fast spatial variations of the original electron fields 
contained in the rapidly oscillating e^^^''^ factors. The 
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low-energy kinetics of the helical edge electrons is then 
encoded by the one-dimensional Dirac Hamiltonian 



-ivF / dx 



(1) 



with vp the Fermi velocity. Away from half filling of 
the ID band of edge states, time-reversal invariance con- 
strains the possible e-e interaction processes to dispersive 
(~ gd) and forward (~ gf) scattering 0, Close to the 
Fermi points zLkp these interactions are given by 

Hd = gd f dx^l^^^l^^, Hf = ^ [ dxViV'a^lV'a, (2) 



respectively, with repeated spin indices a summed over. 
While the forward scattering only leads to a trivial ve- 
locity shift, the dispersive scattering, controlled by Hd, is 
an important feature of a QSH edge In the present 
problem it will be shown to have dramatic effects. 
We now add a Rashba SO interaction [loj . 



Hp, 



dxa{x)¥jx)ayk,^Jx)+R.c., (3) 



where a{x) is a spatially varying coupling, is a Pauli 
matrix, and is the wave number along the edge. 
Rewriting Hn in terms of the slowly varying fields ip^ 
and V'i, with kx~~idx, one obtains 



Hb = dx a{x) ( ((93;'0|)'0^ - ^jldxip^ ) e 



-2ikFX 



-H.c. (4) 



With an eye towards QSH physics in a HgTe QW, a 
few comments are in order. First, the lowest bulk con- 
duction band {Hi) in an inverted HgTe structure is non- 
parabolic, implying that the spatially averaged effective 
bulk Rashba coupling, a* call it, is k dependent [l5[. As 
an estimate of the spatially averaged Rashba coupling 
(a(x)) for the edge states in Q we shall take the size 
of a* in the range corresponding to the wave numbers 
where the two edge bands join the Hi band. Second, one 
should realize that the coupling a{x) depends on several 
distinct features of the QW, most importantly the ap- 
plied gate electric field, the ion distribution in the nearby 
doping layers [l2j , and the presence of random bonds at 
the two QW interfaces As we shall need no details, 
we here treat a{x) as a phenomenological parameter, and 
write it as a{x) = {a{x)) + J2n'^i^n)e^''"^ , with d(fc„) 
the Fourier modes of the zero-mean random contribution 
from the dopant ions and interface bonds, here taken to 
obey Gaussian statistics [l6| . 

With these preliminaries, let us now explore the effect 
of the Rashba coupling on the edge dynamics. As ex- 
pected, the spatial average {a{x)) is seen to leave no trace 
in Hji as the corresponding terms in the integrand oscil- 
late rapidly and average out upon integration. Hence, 
to lowest order, a perturbation with a uniform Rashba 



interaction has no influence on the low-energy dynam- 
ics of the edge states 17[. In fact, the invariance of Hr 
under time reversal T, together with the = — 1 prop- 
erty of a single-electron state, implies that the lowest- 
order effect produced by Hji for any a{x) can at most be 
of 0{a^). For noninteracting electrons, power counting 
on Eq. (|4]) shows that an 0{a^) process is irrelevant [in 
rcnormalization-group (RG) sense], implying robustness 
of the edge states against perturbations with a Rashba 
interaction, even when spatially fluctuating. 

To find out how this picture may change in the pres- 
ence of e-e interactions, we bosonize the Hamiltonian 
H = Ho + Hd+ Hf + Hr, defined by Eqs. ©, @ 
and (HI, and put = t]^ exp {iy/n\6 + 9]) / ^/2JTK and 
-0^ = rj^exp — 0]) /^/2ttk [18[. Here </> and 

6 are dual bosonic fields satisfying dt<j> = vdxd with 
v= [{vF+gf /t^Y — {gd/T^YY^'^ , and 77^^ are Klein factors 
satisfying {ry-i-, 77^} = 0. The microscopic cutoff k is given 
by the penetration depth of the edge states, k sa vp/W , 
with W the bulk band gap (in units with h=l). Absorb- 
ing a factor of i in 77^, one obtains 

Ho + Hd + Hf^v I dxy^{dx(t>f + ^{dxef^ , 

Hn = I dx7]{x){dxe)e''^''' + II.C., (5) 



with r]{x) = J2n a(A:„ — 2fcF)e and K = [{nvp + gf — 
ffrf)/(7rwF +5/ 

Having obtained the theory on bosonizcd form, Eq. 
([5]). we pass to a Lagrangian formalism by a Legendre 
transform of the Hamiltonian. We use that H = \fKdxQ 
serves as conjugate momentum to (I)/\/K and integrate 
out H from the partition function Z to arrive at 



Z 



Vipe 



(6) 



with the Euclidean action 



SW] dxdT(^{dr^f + v{dxVf) 



/ dxdr if (2 

2ttk J ^ ^ 



H.c). (7) 



Here ^{x) = 1/{4Kvk) ^„ „, d(fc„ - 2kF)a{kn' - 2kF) x 
^i{k„+k„,)x^ A/f EE VTe^, and ip = 4)/s/K. Note that by 
integrating out H, terms linear in rj become proportional 
to a total time derivative and hence vanish. Terms pro- 
portional to Tj if contribute only an immaterial constant. 

By averaging over the randomness in S[ip\, using the 
Gaussian statistics P[^(a;)] = exp[— / da; ^* 
so that (^(x)) = and {C{x)i{x')) = D^S{x - x'), the 
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replica method [18| yields the disorder-averaged action 
- 77r^y2 / d-^drdr' cos[XK^ab{x,T,T')] 

[ZTTK) J 

where a, 6=1, n are replica indices with ^abix,T,T') = 
<^a(x,T) - (/>6(a;,r'), and = n,/{8TrK^v^) {Rca^{k)) 
with rii the composite density of the dopant ions and in- 
terface bonds that produce the randomness in the Rashba 
coupling The second-order RG equations of D^,v, 

and K, generated by the scaling (r, x) — ?> (r, cc) exp(— ^) 
{£ > 0), are given by dD^/d£ = (3 - 8K)D^,dv/d£ = 
-2vKD^ and dK/di = -2K^D^. It follows that the 
Rashba coupling grows under renormalization when K < 
Kc = 3/8, driving a transition to an Anderson- type lo- 
calized state. The value Kc = 3/8 was identified in 
Refs. @] as the critical value below which correlated 
backscattering in the presence of quenched disorder may 
cause localization of the helical edge modes. Our anal- 
ysis shows that a Gaussian distributed random Rashba 
coupling yields a precise microscopic realization of this 
type of process. 

Importantly, to find out whether the edge electrons of a 
given experimental sample do become localized, one must 
test for the condition ^[q^, < L, with ^[q^, the localization 
length and L the length of the sample. For this purpose 
we rewrite the replica RG equations for and K on the 
form of Kosterlitz-Thouless (KT) equations, dz\\/dl = 
—z\ and dz±/d£ ~ — 2j|Zj^, with z\\ = 4K — 3/2 and z± = 
SKyi^D^. We renormalize z±{l) until i = £* , defined 
by D^{£*) ~ 1. The renormalized localization length 
6oc(^*) is of the same size as the cutoff k, and since all 
lengths in the system scale with a factor exp(— the 
true localization length of the system is given by ^loc ~ 
Kexp(r) [H. 

To make an estimate of ^loc for a HgTe QW we 
must put numbers on our parameters. Starting with 
the random Rashba coupling arandl^); '^c have that 

\/ ^'-'^rand (*'')) ~ i^i^)) fo'" ^ QW with a zinc-blcnde 
lattice structure [l^. This, together with the estimate 
h{a{x)) w ha* « 5 X 10~^^ eVm at the wave numbers 
at which the edge bands join the Hi band [l^, suggests 
that « 1.5 X 10~^, keeping in mind that a normal 
distributed a{x) yields a normal distributed a{k) with 
(Rea^ik)) = {a^{x)) /2, and that (a^) ^3 (a^)^ . As for 
the unrenormalized interaction parameter K (from now 
on denoted Kq), we use that Ko~{l + X\n{d/s))-^/^ for 
9(1 — Sf ~ViQ~0)/fiVF with V{q) the Fourier transform 
of the Coulomb potential, d being the distance to the 
nearest metallic gate, s = max [ k, dqw ] with dqw the 
thickness of the QW, and A = 2e'^ / {Tr'^eeohvp) with e the 
permittivity of the doping and insulating layers between 
the gate and the QW. By varying the composition and 



the thickness of the layers. Kg can take values from « 1 
(fully screened electrons) down to ^ 0.1 (strongly inter- 
acting limit), with 0.5 < Kq < 0.55 in the experiments 
reported in Ref. i Taking s« 10 nm, n, sa 107m, 
and V vf ~ 5x 10^ m/s 0, we have all the data needed 
for estimating how ^loc depends on Kq. The result is plot- 
ted in Fig. [1] revealing that the edge of a micron-sized 
sample localizes for Kq < 0.25. 
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FIG. 1: (Color online) The edge localization length ^loc for 
different values of the interaction parameter Ko . The dashed 
line marks the length of a micron-sized HgTe QW sample [^. 

When ^loc < L, the edge electrons are localized up to 
temperatures of the order of ^t'/(fcs^ioc)i with an expo- 
nential decrease of the conductivity for lower tempera- 
tures [31. Thus, the edge remains insulating at the low 
temperatures (30mK— 2K) at which experiments are typ- 
ically carried out Q . 

Having explored the Rashba interaction with a spa- 
tially random coupling, one may ask how the QSH edge 
modes respond to a periodically modulated coupling. 
The question becomes particularly interesting in light of 
a proposal by Wang [24I to use a gate-controlled Rashba 
couplin g as a current switch in a quantum wire. As shown 
in Ref. |22| . a periodic modulation of a Rashba interac- 
tion, produced by a sequence of equally spaced nano- 
sized gates, makes the electrons backscatter coherently 
and the current gets blocked. By decharging the gates, 
the current is free to flow again. If this blueprint for a 
spin transistor can be copied over to the edge of a QSH 
insulator, one may envision a device enabling fast and 
efficient control of a dissipationless current. 

As a simple model for a periodically modulated Rashba 
couphng we take a{x) = Acos{Qx) in Eq. ([4]) By 
inspection, the largest effect from this modulation is ob- 
tained by choosing Q = 2kp, as this will produce terms 
in the integrand where all oscillations are carried by the 
slowly varying fields tl'-\,i (with 2fci? corresponding to a 
wavelength of roughly 5—10 nm in a HgTe QW fl^). 
By considering Kq values corresponding to ^loc > L (cf. 
Fig. 1), we can neglect the intrinsic disordering of the 
Rashba coupling. Repeating the procedure that took us 
from Eqs. Ill), (0), and (g]) to the action in Eq. but 
now with r]{x) replaced by the amplitude A, we obtain 

S[f] = ^ 1 dxdT(^^{dripf + v{dxf) - 2g cos{XKf)y 
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with g = / {AttK k'^v). This action is that of the wcU- 
known quantum sine-Gordon model. The RG flows of g 
and K are governed by the KT equations dzu/d£ = —z\ 



and dzx^jdl 



-zi\z±, with 2|| — 4:K — 2 and z_i_ — 



igy/CK^, where C is a K-dependent constant |18| . When 
Kq < 1/2, the mass term in the sinc-Gordon action is 
seen to be RG relevant. In this regime backscattering 
generates a dynamical mass for the edge modes, turning 
the edge to a Mott insulator. Whether the correspond- 
ing mass gap A, measured from the Fermi level, is large 
enough to block the current depends on its size relative to 
the available thermal energy. To make an estimate for the 
case of a HgTc QW of length 1 /im [20 fim], we calculate A 
for three different values of the Rashba amplitude, tuning 
up the interaction parameter Kq until the renormaliza- 
tion length Kexp(^*) reaches 1 ^m [20 ^m], with exp(^*) 
the scale factor for which backscattering processes start 
to dominate (see Fig. [5]). When Kexp(^*) « 1 /.im, 
A « 4.5 meV for all three amplitudes, corresponding to 
a temperature of around 50K. 

6 
4 



-hA = 7 X 10""cVm 
-hA = 5 X 10""cVm 
-7iA = 3 X 10-"cVm 
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FIG. 2: (Color online) The gap A for different Rashba am- 
plitudes hA and values of Ko The circles and squares mark 
the smallest gaps for HgTe QW samples of length 1 /im and 
20 /im, respectively. 

The thermally activated conductance is given by G = 
{2e^ /h)e^'^/'^^'^ . To prevent thermal leakage of current 
through the proposed switch, one may thus need to cool a 
micron-sized sample to temperatures in the range 10 — 20 
K. For these temperatures our result provides a proof-of- 
concept of using an " on-off" modulated Rashba coupling 
as a switch for QSH edge currents. It should be interest- 
ing to repeat our analysis for QWs based on the ternary 
Heusler compounds |24l |. once these have been pu t to- 
gether and characterized. As discussed in Ref. |2J, the 
diversity of Heusler materials opens wide possibilities for 
tuning the bulk band gap and setting the desired band 
inversion. This may allow also for enhanced Rashba in- 
teractions at the edge. The prospect for a QSH-based 
spin transistor makes this a promising path to explore. 

In summary, we have carried out an analysis of the 
combined effect of Rashba spin-orbit and e-e interactions 
on the edge dynamics in a QSH state. The spatial dis- 
ordering of the Rashba coupling, intrinsic to a quantum 
well structure that supports a QSH state, is found to lo- 
calize the edge electrons for sufficiently strong e-c inter- 
actions. Likewise, a periodic modulation of the coupling 



— realizable by placing a configuration of equally spaced 
nanosized gates on top of the quantum well — localizes 
the electrons already at intermediate strengths of the e-e 
interaction. While a practical implementation may have 
to await advances in nanoscale technology, our result sug- 
gests that a Rashba-controlled switch for the dissipation- 
less edge current in a QSH device is a viable possibility. 
On the more fundamental side, trying to characterize the 
transition from the Mott insulating phase caused by a 
periodic Rashba modulation to the Anderson-type local- 
ization due to the interplay of e-e interactions and the 
disordered Rashba coupling remains a challenge. 
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